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[ Within the kinetic approach velocity and charge distributions of grains in stationary dusty plasmas 
are calculated and the relations between the effective temperatures of such distributions and plasma 
parameters are established. It is found that the effective temperature which determines the velocity 
grain distribution could be anomalously large due to the action of accelerating ionic bombarding force. 
Ch The possibility to apply the results obtained to the explanation of the increasing grain temperature 
in the course of the Coulomb-crystal melting by reduction of the gas pressure is discussed. 
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^ ■ Recently much attention has been payed to theoretical studies of various problems of dusty 
, plasma physics associated with grain dynamics and grain charging (formation and melting of 

\ dusty crystals, influence of charging on effective grain interaction, dust-acoustic wave excitation, 
] effect of grain charging on fluctuations and electromagnetic wave scattering in dusty plasmas, 
00 ' etc.). In such studies it is convenient to treat the grain charge as a new variable (as was done 
for the first time in Ref. [1]). This makes it possible to statistically describe the grain charge 
distribution on equal footing with the spatial and velocity grain distributions. Obviously, it 
is very important to know what are the stationary (quasiequilibrium) grain distributions and 
Q\ , what is the relation of these distributions to plasma parameters. In spite of the fact that 
i statistical descriptions of dusty plasmas have been already used in many papers, as far as 

■ the authors of this letter know neither grain charge, nor velocity distributions for grains were 
studied within a consistent kinetic approach. Usually, the problem is avoided by neglecting 

' the thermal dispersion of grain velocity and charge. In many cases this is a rather reasonable 
approximation, but it could not be valid when the properties of the grain subsystem and its 
^ dynamics are concerned. 
X ■ The purpose of the present paper is to describe stationary velocity and charge distributions 

■ of grains in dusty plasmas in the case of grain charging by plasma currents and to determine 
the dependences of effective temperatures on plasma parameters. We study dusty plasma 
consisting of electrons, ions, neutral molecules and monodispersed dust particles (grains) as- 
suming that every grain absorbs all encountered electrons and ions. Such collisions we define 
as charging collisions. Collisions in which plasma particles do not touch the grain surface we 
call Coulomb elastic collisions . Notice that the cross-sections of charging collisions are also 
determined by the Coulomb forces along with the geometrical size of grains. 

Using the microscopic equations for dusty plasmas and the relevant BBGKY-hierarchy [2] 
it is possible to show that in the case of dominant influence of charging collisions the kinetic 
equation for the grain distribution function fg{X, t) = /g(r, v, q, t) {q is the charge of the grain) 
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where (Tpo.(g, v) is the cross-section for charging: 
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6{x) is the Heaviside step function, a is the grain radius, fa{r,v,t) is the plasma particle 
distribution function normalized by the particle density n^, Sv^ = (^cr/^(/)v' is the grain 
velocity change due to the collision with a plasma particle, subscript a labels plasma particle 
species, the rest of the notations is traditional. Eq. (1) could be introduced also on the basis of 
physical arguments as was done in Refs. [3,4]. In fact, the right-hand part of Eq. (1) describes 
the balance between the grains outcoming from the phase volume element and those incoming 
to the same element due to charging collisions. 

Taking into account the smallness of and Sv^ it is possible to expand the right-hand part 
of Eq. (1) into a power series of these quantities. With the accuracy up to the second order 
Eq. (1) in the stationary isotropic and homogeneous case is reduced to 
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where D\\, /3, Q, 7 and / are the Fokker-Planck kinetic coefficients generated by charging 
collisions and given by 
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The quantities D\\{q^v) and Q{q,\) characterize the grain diffusion in the velocity and charge 
space, respectively, P{q,v) and 7(9', v) are the friction coefficients which determine the bom- 
bardment force Fb(g, v) = —mg(3{q,v)v associated with charging collisions and the correction 
to this force (5Fb(g',v) = — mg7(g, v)v due to the mutual influence of the charge and velocity 
grain distributions, I is the grain charging current. Deriving the relation for f3{q,\) we omit 
the terms of higher order in {uifj/mg) associated with the tensor nature of the diffusion coef- 
ficient in velocity space (contribution of the transverse diffusion coefficient). With regard for 
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the fact that v)| — > oo at Cg. — > and v)/D||(g, v)| — > oo at {m„/mg) —>■ 0, it 

is possible to show that the asymptotical solution of Eq. (3) can be written as 

where 
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6q{v) = q-q(v), 

q{v) is the stationary charge of the grain moving with the velocity v, given by the equation 

l{q{v),v)^0, (7) 

Z is a normalization constant, £ and A are small functions. Substitution of Eq. (5) into Eq. (4) 
leads to 
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Eqs. (5)-(8) give the asymptotically exact solution of Eq. (3) at (mgecr/mcrq) 00. Further 
estimates require the explicit form of the kinetic coefficients. Assuming that plasma parti- 
cle distributions are Maxwellian, one obtains the following stationary grain distribution with 
accuracy up to the zeroth order in {qmi/eeiTig): 
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UQg is the averaged number density of grains. Here, we use the notation 



The quantity go is the equihbrium grain charge of stationary particles satisfying the equation 
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For typical values of plasma parameters in dusty plasma experiments {t + z > \) and Zi = 1 
we have: 

In such case the thermal variation of the grain charge |g — gop is of the order of aT^ and 



|g — go 1-2 ~ %\l ^1/ dT^. This means that at weak plasma coupling defined with the grain size 
{el/aTe <S 1) the effective temperature of the grain thermal motion reff(g) = Tes reduces to 

Teff ^ 2T,^ (14) 

and 



D{q, v) ~ Do (l + ^) , g(g, v) ~ Qo{t + z){l + Z,). 
Thus, in such case 
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This distribution describes the equilibrium Maxwellian velocity distribution and the Gibbs 
grain charge distribution with the temperatures Tcs and Tcs respectively. In fact, the electric 
energy of charge variations of the electric capacity a is equal to (g — go)^/2a and thus, the 
charge distribution described by Eq. (15) can be interpreted as an equilibrium distribution with 
effective temperature Tcff-. At t < 1, 2; < 1 the effective Teg exceeds the electron temperature. 
The resulting velocity distribution is described by the effective temperature Tcff. Even in the 
case of neutral grains {z = 0) this temperature is equal to 2Tj. The presence of the factor 2 is 
associated with plasma particle absorption by grains. 

Charging collisions are inelastic and a part of the kinetic energy of the ions is transformed 
into additional kinetic energy of the grains. This is the difference between the case under 
consideration and conventional Brownian motion where the velocity distribution is described by 
the temperature of the bombarding light particles. Eq. (14) shows that the effective temperature 
of thermal grain motion could be anomalously high at z ^ t. Physically it can be explained 
by the decrease of the friction coefficient with increase of grain charge 

4 V2^g)A.S.(l-g^/3o(l-i 

The reason is that the difference between the fluxes of ions bombarding the grain surface 
antiparallel to the grain motion and parallel decreases with the charge increase due to the 
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specific properties of the ionic charging cross-section, which charge- dependent part is larger for 
ions moving with smaller relative velocities (i.e. in parallel direction). The condition z — t 
corresponds to the zero value of the friction force. 

Eq. (3) and its solutions (5), (9), (15) were obtained under the assumption that the Coulomb 
elastic collisions could be neglected. In order to take elastic collisions into consideration Eqs. (1), 
(3) should be supplemented by the appropriate collision terms, for example, by the Landau, or 
Balescu-Lenard collision integrals. We use the Balescu-Lenard collision integral in the Fokker- 
Planck form which in the case under consideration (isotropic spatially homogeneous stationary 
distribution) can be written as 
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^iD\\c{q,^^)fg{v,q)) + vPc{q,yr)f,{q,v) , (16) 

where D\\c{q,v) and Pc{q,^) are the Fokker-Planck coefficients related to Coulomb elastic 
coUisions (see, for example, [5], Chapter 8). With the accuracy up to the dominant logarithmic 
terms (in this approximation Eq. (16) is reduced to the Landau collision term) such coefficients 
can be reduced to 
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In Eqs. (16), (17) we again neglect the contribution of the transverse part of the diffusion coeffi- 
cient which gives a correction to /3c(?, v) of higher order in {rricr/ nig) and we disregard the grain- 
grain Coulomb collisions, assuming the grain density to be small {ug < ni{Zi/ ZgY{Sg/ SiY/'^{Tg/Ti)). 
We introduced also the Coulomb logarithms In A^- for each particle species. Usually these quan- 
tities are estimated as InA^. = ln(A;inax/^D)! where k£, = r]j" = {Y^i'^'^^ai^'a/TaY^'^ and /Cmax is 
the inverse distance of closest approach between colliding particles, 

m„v^ 3T„ _i 

kraa^a ~ T ~ 1 T = ^La (18) 

Kq\ Kql 

(rLo- is Landau length). However, in the case of plasma particle collisions with finite-size 
grains this estimate could be invalid, since at Tlo- < a the Coulomb logarithm will include the 
contribution of collisions with particles reaching the grain surface, i.e. charging collisions. 

An approximate modification of A^. is achieved by treating InA^ as a logarithmic factor 
appearing in the momentum transfer cross-section for Coulomb collisions. In the case of finite 
size grains one obtains the following logarithmic factor 



m Ao- = m sm / sm 
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where Xmaxa and Xmina are the scattering angles related to the minimal and maximal impact 

parameters ftmino- and ftmaxo- by the Rutherford formula. Obviously, femino- should be determined 
from the condition that the distance of closest approach is equal to a implying 
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bmma = a\l —e{l — ). (19) 
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Concerning the quantity ^maxo-, it is reasonable to put 6maxcr = td + o instead of 6maxcr = ^"d, 
since in the case of a finite size grain its screened potential is given by the DLVO-potential 

$(r) = ^(l + A)-ie-('-«)A^, 

rather than the Debye potential. 
As a result we have 



In A,- = - In 



2 {ru + a) 
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In Ae = - < 
^ 2 



As is seen, at ru ^ ''c the ionic Coulomb logarithm can be a small quantity in contrast to the 
case of ideal plasmas. 

Comparing Eqs. (3) and (16) it is easy to see that in order to take elastic Coulomb collisions 
into account it is sufficient to make the following replacements in the obtained solutions 

D\\ {q, v) ^ D^i {q, v) = L>|| (g, v) + L>||c(g, v) 
P{q, v) ^ P{q, v) = P{q, v) + Pc{q, v). (21) 

In the case of weak plasma coupling {el/aTg <^ 1) 
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(3{q,v)^(3o[l-- + 2-^\nA,j, (22) 

Thus, the correction produced by the elastic collisions could be of the same order as that due 
to charging collisions. The condition for dominant influence of charging collisions is 
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which can be realized at small values of z/t, or at z\t ^ rjj/a^lru S> r|)/a). 

Rigorously speaking Eq. (16) and thus Eq. (22) are definitely valid in the case of weak 
coupling plasmas {ru <C r^) since this is the condition of the derivation of the Balescu-Lenard 
(or. Landau) collision term. However, it is possible to expect that actually the domain of 
vahdity of Eqs. (16), (22) is not too strongly restricted by such condition. This assumption is 
in agreement with the direct calculations of the friction coefficient (Coulomb collision frequency) 
in terms of the binary collision cross-sections. Beside that, as it was shown in Ref. [6], in the 
case of strong grain-plasma coupling the infiuence of the Coulomb collision is also small and the 
kinetic equation is reduced again to Vlasov equation. This means that fluctuation evolution 
equations, which solutions determine the explicit form of the Balescu-Lenard collision term, are 
the same as in the case of weakly coupled plasmas and thus Eq. (16) continues to be valid. 
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The new kinetic coefficients give the foUowing effective temperature for thermal grain motion 

2(i + | + |ilnA,) 

l-f + 2f^lnA, ^ ' 

i.e. elastic collisions can produce a saturation of the grain temperature. However, in the case 
of dominant influence of charging collision Tpfj can be still anomalously large. This fact can be 
used for a qualitative explanation of the experimentally observed grain temperatures which are 
usually much higher than the ion temperature, Tg ^ Ti (see, for example [7,8], Tj ~ 0.1 eV, 
~ 4-=-40 eV). Finally, we point out that the obtained results can be modified also for the case 
of a plasma with a neutral component. It is possible to introduce an additional collision term 
along with the term (16). Since the collision integral describing elastic collisions of neutrals 
with grains also can be represented in the Fokker-Planck form (it follows from the Boltzmann 
collision integral) the presence of neutrals results in new additions to D\\ and /3, namely 

i5.(..)^..(.n4in^.^(:::^)-'-^(l)T 

^ A (l - i . Sfjln A. + 2^ (^)"' (1)"^) . (24) 

As a result the effective temperature is modified into 
(l + f + glnA. + ^fe)^^^(t)^^^) 
(l-f + 2#lnA. + 2H^(^)^ (t) ) 

According to Eq. (25) the effective temperature increases with decreasing neutral density. The 
infiuence of neutral density changes on the effective temperature would be especially important 

2 

at 1 - f + 2f7 In Ai ^ 0. In such decrease of the neutral gas pressure can produce an 

anomalous growth of T^s- That is in qualitative agreement with the experimental observation 
of melting of dusty crystals by reduction of the gas pressure [7,8]. 

The obtained results show that stationary velocity and charge grain distributions are de- 
scribed by effective temperatures different from those of the plasma subsystem. These effective 
temperatures are determined by the competitive mechanics of collisions: grain-neutral colli- 
sions and clastic Coulomb collisions result in the equalization of the effective temperature to 
the temperature of neutrals, or ions, respectively, while charging collisions can produce anoma- 
lous temperature growth. That could be one of the main mechanisms of grain heating. 
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